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Abstract 


The  vortex-lattice  method  is  employed  to  calculate  the  inviscid 
flow  and  the  unsteady  development  of  separated  vortex  sheets  over  a  pro¬ 
late  spheroid.  An  approximate  boundary-layer  method  based  on  the 
assumption  of  local  similarity  is  used  to  calculate  the  line  of  open 
separation.  The  two  schemes  interact  through  the  line  of  separation, 
which  is  allowed  to  displace  as  the  wake  grows.  The  shape  of  the  wake 
is  calculated,  and  pressure  distributions  are  compared  with  experimental 
data.  The  interaction  scheme  predicts  very  accurately  the  pressure  on 
the  attached  flow  region  but  displays  unexpected  special  variations  over 
the  separated  region. 


1.  Introduction 


Aerodynamic  problems  involving  massive  regions  of  separated  flow  in 
three  dimensions  can  be  treated  today  numerically  only  by  methods  of 
discrete  vortex  dynamics.  Navier-Stokes  codes  would  require  incredible 
amounts  of  computer  time  and  space  to  generate  the  solution  to  a  practi¬ 
cal  problem.  In  two  dimensions,  a  variety  of  problems  have  been  con¬ 
sidered  using  discrete  vortex  methods;  for  example,  the  stability  of 
free-shear  layers  and  the  wake  structures  of  thin  plates  or  cornered 
bodies.  More  related  to  the  present  contribution  are  works  considering 
the  problems  of  separation  over  bluff  bodies;  the  most  popular  case 
being  the  circular  cylinder. In  fact,  the  first  efforts  to  consider 
separated  regions  over  elongated  three-dimensional  bodies  were  based  on 
the  cross-flow-plane  analogy  which  employs  two-dimensional  solutions  in 
cross  sections  of  the  body.  Two-dimensional  methods  have  been  refined 
to  include  redistribution  of  vorticity,  models  for  viscous  dissipation, 
etc. 

In  three-dimensional  problems,  one  set  of  such  methods  is  known  as 
"vortex-lattice  methods."  Problems  involving  sharp-edge  separations 
have  been  treated  in  this  way  with  considerable  success  for  steady®'^ 
and  unsteady^*^"^^  flows. 

Three-dimensional  flows  involving  separation  over  smooth  surfaces 
have  been  considered  only  very  recently.  Fiddes^^  expanded  on  a  vis- 
cous-inviscid  interacting  method  to  solve  the  problem  of  three-dimen¬ 
sional  separation  over  a  cone.  This  is  an  elegant  mathematical  thesis 
and  a  revealing  contribution  limited  to  slender  body  configurations  and 


Thrasher^^  studied  the  flow  over  a  semi-infinite  body  with  a  tan¬ 
gent-ogive  nose  and  cylindrical  afterbody.  He  employed  an  iteration 
version  of  a  vortex-lattice  scheme  to  solve  the  problem  and  assumed  that 
the  separation  lines  coincided  with  generators  of  the  body.  Almosnino 
and  Rom^^  have  considered  the  flow  about  a  very  similar  semi-infinite 
body  and  employed  again  an  iteration  method  to  achieve  convergence  of 
the  wake  vortices.  However,  they  modeled  the  separated  vortex  sheets  by 
only  four  vortex  lines  on  each  side  emanating  from  very  narrow  segments 
aligned  again  with  generators  of  the  body. 

In  all  invesitgations  of  three-dimensional  bluff-body  separated 
flows,  the  body  was  assumed  to  be  semi-infinite.  Moreover,  with  the 
exception  of  the  degenerate  case  considered  by  Fiddes,^^  all  other 
investigators  assume  that  the  position  of  separation  is  known  and  that, 
in  fact,  it  coincides  with  a  meridional  line  of  axisymmetric  bodies.  In 
the  present  paper,  we  report  tn  our  initial  efforts  to  relax  the  above 
restrictions.  To  this  end,  we  have  coupled  the  potential  flow  calcula¬ 
tions  with  an  approximate  method  of  calculating  the  three-dimensional 
boundary  layers.  Moreover,  we  considered  the  unsteady  development  of 
the  wake  and  allowed  the  viscous  and  inviscid  fields  to  interact.  The 
line  of  separation  over  a  finite  body  was  readjusted  at  each  time  step 
and  was  found  to  tend  to  the  position  predicted  experimentally  for  a 
fully  developed  flow. 

Three-dimensional  laminar  boundary  layers  have  been  calculated  by  a 
variety  of  methods  in  the  past  decade.  We  reference  here  the  most 
recent  contributions  of  Cebeci  et  al.,^^  Patel  and  Baek,^^  Tai,^®  and 
Ragab.^^  All  four  papers  contain  extensive  reference  to  earlier  contri¬ 
butions.  The  present  authors  have  developed  an  aporoximate  method^*^ 


which  is  described  briefly  here  and  is  employed  for  ,he  calculation  of 
the  separation  line. 


» 


A  significant  feature  of  the  flow  and  a  controversial  point  is  the 
development  of  open  separation  lines.  This  phenomenon  was  described  by 
Maskell^^  and  very  recently  discussed  by  Tobak  and  Peake^^  and  Wang.^^ 
In  the  present  report  ve  reference  briefly  the  evidence  of  our  experi¬ 


mental  work  in  favor  of  open  separation^^  and  employ  our  analytical 


method  which  is  in  agreement  with  this  idea. 


2.  Potential  Flow  Calculations 


The  outer  flow  and  the  separated  vortex  sheet  are  predicted  by  the 
vortex-lattice  method.  This  is  based  on  6iot-Savart ' s  law  which  stipu¬ 
lates  that  the  velocity  vector  V  at  position  r  induced  by  a  vortex  fila¬ 
ment  of  strength  r  is  given  by 


nr)  -  ^  1 


(1) 


r-s 


V(r)  =  e  ^  (cose^-cose^) 


where  the  variable  s  is  a  position  vector  of  points  along  the  vortex 
line.  An  equivalent  formula  for  the  present  application  gives  the  velo¬ 
city  induced  by  a  vortex  segment  of  finite  length.  With  quantities 
defined  in  Fig.  1,  the  velocity  is  expressed  as 

(2) 

where  e  is  a  unit  vector  normal  to  the  plane  defined  by  the  point  of 
interest  and  the  vortex  segment. 

The  velocity  field  calculated  by  the  superposition  of  quantities  of 
this  form  is  irrotational  except  along  the  vortex  lines  which  represent 
singularities  of  the  solution.  Appropriate  distributions  of  vorticity 
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along  fixed  or  moving  grids  can  be  used  to  generate  solutions  of  incom¬ 
pressible  flows  in  two  or  three  dimensions  with  or  without  separation. 

As  a  test  for  three-dimensional  bluff  bodies,  we  considered  three  basic 
problems  which  accept  exact  closed-form  analytical  solutions,  readily 
available  for  comparison:  the  flow  about  (a)  a  sphere,  (b)  a  prolate 
spheroid,  and  (c)  a  general  ellipsoid.  To  solve  these  problems,  vortex- 
lattices  have  been  defined  with  their  nodes  on  tne  surface  of  the  body 
as  shown  in  Figs.  2  and  3.  The  strength  of  each  vorticity  segment, 
usually  referred  to  as  branch  circulation  is  denoced  by  the  symbol  :. 
Most  investigators  working  in  this  area^’^* express  Eg.  (2)  in 
terms  of  the  loop  circulation,  G,  around  a  panel  of  the  vortex  lattice. 
The  branch  circulation  :  is  simply  the  difference  of  the  loop  circula¬ 
tions  of  its  two  neighboring  panels,  Gj^  and  02* 

r  =  G^-G^  (3) 

The  branch  circulations  thus  defined  automatically  generate  a  diver¬ 
gence-free  vortex  sheet,  which  is  a  classical  property  of  vorticity. 

For  convenience  of  the  notation,  consider  a  case  of  N  panels  with  N 
unknown  loop  circulations  Gj,  j  =  1,...N.  These  quantities  can  be 
specified  by  imposing  the  no-penetrat ion  condition  at  N  points  on  the 
body.  We  chose  here  the  "average  control  points,"  namely,  the  centroids 
of  the  panels.  However,  we  project  these  points  on  the  actual  surface 
of  the  body  under  consideration.  The  influence  matrix  is  first  con¬ 
structed.  An  element  of  this  matrix  is  the  value  of  the  normal  velocity 
at  tne  ith  control  point  due  to  the  vorticity  distributed  around  the  jth 
panel  with  a  un'^t  loop  circulation.  The  no-penetration  condition  then 
is  simply  expressed  as 


where  is  the  freestream  and  n^-  is  the  normal  to  the  body  at  the  point 
i.  The  flowfield  about  the  body  under  consideration  can  be  generated  by 
solving  the  system  of  equations  for  the  quantities  Gj. 

Typical  results  obtained  by  such  calculations  are  shown  in  Figs,  d 
and  5  where  ^e  compare  the  velocity  distribution  along  a  meridional  cal¬ 
culated  by  the  exact  method  and  the  vortex  lattice  method.  The  results 
displayed  in  these  figures  correspond  to  an  optimum  choice  of  panel 
dimensions.  Surprisingly,  the  accuracy  of  the  calculation  is  improved 
by  adding  elements  in  the  axial  direction  alone.  Excellent  results  were 
obtained,  for  example,  with  27  panels  in  the  axial  direction  and  only  6 
or  even  just  4  in  the  circumferential  direction. 

Separated  flows  can  be  calculated  by  approximating  the  separated 
vortex  sheets  again  by  vortex  lattices. This  is  discussed  in 
detail  in  Section  4.  In  the  present  case  we  chose  to  consider  impul¬ 
sively  started  flows.  The  steady  state  will  then  be  obtained  as  a  limit 
for  large  times.  At  present,  and  mostly  due  to  limitations  in  computing 
time,  we  have  not  obtained  large-time  solutions.  However,  the  present 
results  display  some  very  encouraging  characteristics. 

The  Kel vin-Helmholtz  theorem  dictates  that  vorticity  be  transported 
with  the  local  velocity.  This  condition  automatically  generates  a 
force-free  wake  and  is  used  here  in  a  manner  similar  to  the  one  des¬ 
cribed  in  earlier  references. The  main  difference  here  is  that  the 
wake  is  initiated  along  a  line  of  separation  which  is  determined  by  a 
viscous  calculation  and  is  free  to  readjust  at  each  time  step. 

The  shed  vorticity  in  the  form  of  two  rows  of  panels  at  the  root  of 
the  separated  vortex  sheets  are  determined  at  each  time  step  from  vis- 
cous-inviscid  interaction  as  described  in  Section  4.  In  the  next  time 


step  this  row  of  panels  is  convected  away  from  the  body  and  a  new  row  of 
panels  is  generated  in  the  same  way.  The  free  vortex  panels  must  be 
accounted  for  in  the  calculation  of  the  velocity  field.  Equation  (4) 
must  be  modified  to  include  V,  ,  the  velocity  induced  by  the  wake 


y  A.  .G,  =  (V  +  V  )-n. 


(5) 


It  should  be  noted  that  the  matrix  A^-j  is  universal  for  a  certain  body 
configuration,  but  the  calculation  of  the  velocity  although  tedius 
and  time  consuming,  must  be  repeated  at  each  time  step.  Finally,  the 
system  of  Eqs.  (5)  is  nonsingular  and  can  be  solved  with  the  same 
routine  used  to  solve  Eqs.  (4). 


3.  Boundary-Layer  Calculations 

The  boundary  layer  was  calculated  by  an  approximate  method 
developed  earlier  by  the  authors^^  and  repeated  here  for  completeness. 
The  method  is  based  on  the  assumption  that  the  boundary-layer  flow  o'^er 
a  general  body  is  locally  similar  to  the  flow  over  a  special  body  whose 
solution  is  readily  available.  Such  bodies  are  cylinders  of  arbitrary 
cross  section  at  yaw  with  respect  to  the  oncoming  flow. 

Consider  the  flow  about  a  cylinder  in  yaw.  Let  the  x,  y,  and  z 
axes  be  parallel  to  the  generators,  along  the  contour  of  the  body  and 
perpendicular  to  the  generators,  and  normal  to  the  body,  respectively 
(see  Fig.  6a).  Let  u,  v,  and  w  be  the  corresponding  velocity  compon¬ 
ents,  and  U  ,  V  the  components  of  the  undisturbed  freestream.  The 

33  X 

governing  equations  in  their  dimensionless  stretched  form  then  read 


•  •  •  >  •  4  * 


=  0 

3y  3Z 


V  l!i  +  w  ^  = 


2 

3V  .  3V  ,,  dV  3  V 
3y  3Z  dy  ^^2 


(6) 


(7) 


(8) 


Notice  that  there  is  no  pressure  gradient  term  in  Eq.  (7).  Moreover, 
Eqs.  (6)  and  (8)  can  be  solved  independently  of  Eq.  (7).  The  familiar 
boundary  conditions  are 


u=v=w=0  atz=0 


u-U.v-V  asz 


(9) 


(10) 


For  an  outer  flow  given  by 


U  =  U 


V  =  C^y  +  Cjy 

the  solution  was  obtained  in  the  form  of  two-term  expansions 


u  =  U^(go(z)  +  ^  y^g2(z)l 


(11) 

(12) 


V  =  C^y  f^(z)  +  C,y  f,(z) 


3J-  '3^ 


(13) 

(14) 
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and  the  functions  f,,  fj,  g^  and  g-  were  tabulated. 

According  to  our  method,  the  rules  of  correspondence  between  points 
on  the  body  of  interest  and  a  Sears  body  required  matching  of  the 
fol  lowing: 


*•'»  *•  * 


f  f.'f 


a. 


the  reduced  distances  x  along  inviscid  streamlines  scaled  with 
the  total  length  of  each  streamline, 
the  U  velocity  component. 


b. 

c.  the  V  velocity  component,  and 

d.  the  velocity  gradient  along  the  inviscid  streamline. 

For  a  coordinate  system  x,ii  defined  on  a  prolate  spheroid  as  shown 
in  Fig.  6b,  the  four  conditions  read 


X  =  X 

s 

(15) 

3 

v(x,4))  — 

(16) 

u(x,4;)  =  V 

(17) 

v^.l.  =  c  +  3C  x^ 

3S  "*^3  S 

(18) 

where  is  the  reduced  distance  along  the  Sears  body. 

Equations  (15)-(18)  are  solved  at  each  point  on  the  body  of 
interest  for  the  quantities  x^,  Cj^,  C3,  and  V.  They  thus  define  a  Sears 
body  and  determine  a  point  on  this  body.  The  solution  at  this  point  is 
then  given  by  Eqs.  (13)  and  (14)  and  is  assumed  to  hold  at  the  corres¬ 
ponding  point  on  the  body  of  interest. 

The  direction  of  the  skin  friction  on  the  surface  of  the  body  was 
calculated  in  this  way  at  the  control  points  of  the  vortex  lattice 
method.  A  computer  program  was  prepared  for  the  calculation  of  skin 
firction  lines  and  boundary-layer  edge  streamlines.  Integration  was 
started  at  points  equally  spaced  on  a  meridional  line  and  was  marched  in 
both  the  upsteam  and  downstream  direction.  In  this  way  difficulties 
encountered  in  the  neighborhood  of  stagnation  points  were  avoided.  In 
fact,  the  method  provides  a  clear  indication  of  accuracy  at  stagnation 


points  where  it  was  demonstrated  that  ail  lines  pass  within  a  distance 
of  a/1000  of  each  other. 

4.  Viscous-Inviscid  Interaction 

In  the  present  problem  we  allow  the  boundary  layer  to  interact  with 
the  outer  flow  only  through  separation.  The  stength  of  the  vorticity 
shed  at  separation  and  the  location  of  the  line  of  separation  control 
the  development  of  the  free  vortex  sheets  which  roll  up  forming  the 
wake.  The  structure  of  the  flow  at  separation  and  the  mathematical 
modeling  that  allows  the  two  solutions,  the  viscous  and  the  inviscid  to 
interact  with  each  other  is  of  crucial  importance. 

The  topology  of  skin  friction  lines  has  been  studied  extensively  in 
three-dimensional  boundary  layer  flows.  Experimental  evidence  and 
numerical  information  indicate  that  the  skin  friction  lines  merge 
together  along  the  line  of  separation.  Recent  careful  boundary-layer 
calculations^®  point  to  the  direction  that  the  separation  line  is  an 
evelope  of  skin  friction  lines  and  therefore,  within  the  framework  of 
uninteracted  boundary- layer  theory,  a  singular  behavior  should  be 
expected  there.  Thus,  it  is  justified  to  use  the  Sears  model  for  the 
construction  of  local-similarity  solutions,  since  it  is  also  singular 
along  the  line  of  separation. 

It  is  well  known  that  the  vortex  lines  on  the  surface  of  the  body 
are  orthogonal  to  the  skin  friction  lines.  This  implies  that  the 
surface  vortex  lines  meet  the  line  of  separation  at  an  angle  of  90°  and 
lift  off  the  body  surface  also  perpendicular  to  the  line  of  separation. 
Incidentally,  this  is  contrary  to  the  common  concept  of  two-dimensional 
separation  where,  in  principle,  surface  vortex  lines  are  parallel  to  the 
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line  of  separation.  The  authors  and  their  associates  have  obtained  some 
experimental  evidence  that  in  two-dimensional  flows,  wel 1 -organized 
cells  of  cross  flow  are  present  in  the  neighborhood  of  separation.  It 
is  therefore  possible  that  all  separation  lines  are  lines  of  zero  vor- 
ticity.  The  picture  described  above  is  actually  the  limit  of  the 
behavior  of  streamlines  and  vortex  lines  as  the  distance  from  the  wall 
tends  to  zero.  Actually,  the  total  amount  of  vorticity  parallel  to  the 
lines  of  separation  is  not  at  all  zero  at  separation.  In  fact,  it  is 
the  dominant  portion  of  the  vorticity  that  the  boundary  layer  sheds  in 
the  wake. 

A  straightforward  integration  of  vorticity  across  the  boundary 
layer  indicates  that,  within  the  boundary-layer  approximation,  the 
magnitude  of  the  vorticity  flux  at  a  station  of  the  boundary  layer  is 
proportional  to  the  square  of  and  its  direction  is  normal  to  the 
direction  of  the  edge  velocity.  This  fact  has  been  employed  extensively 
in  investigations  of  two-dimensional  flows  and  was  lucidly  explained 
again  by  Sears. The  idea  is  that  in  approximating  the  shed  vorticity 
by  discrete  potential  vortices,  as  shown  schematically  in  Fig.  7a,  the 
strength  of  the  nascent  vortices  must  be  compatible  with  the  vorticity 
contained  in  the  boundary  layer.  For  a  straightforward  derivation  of 
the  classical  formula,  consider  the  circulation  about  the  loop  ACOO' 
shown  in  Fig.  7b. 

r  =  ;  V-dl  =  If  p-dA  (19) 

For  steady  flow,  the  time  derivative  of  the  integral  on  the  right-hand 
side  vanishes,  except  for  a  contribution  along  the  side  00',  which 
allows  vorticity  to  be  convected  outside  of  this  domain. 
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dr 

hF  =  j  uKdy 


And  within  the  boundary-layer  approximation,  a  ~  au/ay  and,  hence, 


dt  2 


where  Ug  is  the  velocity  at  the  edge  of  the  boundary  layer.  This 
condition  is  extended  below  for  three-dimensional  flows  and  is  employed 
to  calculate  the  strength  of  the  nascent  vortices  at  each  time  step. 

Consider  a  three-dimensional  boundary  layer  over  a  solid  surface 
aligned  with  the  coordinate  plane  xy.  In  the  notation  of  Fig.  8  then 
the  flux  of  vorticity  can  be  represented  by  two  components  which  can  be 
obtained  by  similar  integrations  across  the  boundary  layer 
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What  is  of  great  importance  and  interest  here  is  that  nascent 
vorticity  is  parallel  to  the  wall  but  not  necessarily  parallel  to  the 
line  of  separation.  This  is  because  by  definition  the  vorticity  vector 
in  an  attached  boundary  layer  is  parallel  to  the  wall.  The  vorticity 
vector,  therefore,  is  normal  to  the  edge  velocity  which  in  turn  is  not 
normal  to  the  line  of  separation.  If  we  need  to  discretize  the  nascent 
vorticity,  then  we  can  simply  resolve  it  in  components  parallel  and 
perpendicular  to  the  line  of  separation  but  always  parallel  to  the 
surface  of  the  body,  as  shown  schematically  in  Fig.  8.  Let  the  plane  yz 
in  this  figure  be  a  plane  of  separation.  Since  vorticity  is  convected 
with  the  flow,  the  vorticity  vector  is  turned  to  align  itself  with 
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the  flow  direction,  as  shown  again  in  Fig.  8.  This  condition  was 
implemented  in  the  numerical  scheme  employed  here. 

A  discretization  scheme  of  the  free  vortex  sheet  originates  with  a 
ribbon  with  rectangular  elements.  The  surface  of  the  ribbon  was  assumed 
to  be  parallel  to  the  surface  of  the  body  at  an  arbitrary  distance  of 
0.2R  where  R  is  radius  of  the  cross  section  of  the  body  at  the  station 
under  consideration.  More  properly,  this  distance  should  be  a  length 
proportional  to  the  displacement  thickness  of  the  boundary  layer.  The 
outer  nodes  of  the  ribbon  were  allowed  to  be  convected  and  a  new  ribbon 
of  vorticity  was  generated  according  to  the  above  scheme.  In  this  way  a 
free  vortex  sheet  was  generated. 

5.  Results  &  Discussion 

Two  types  of  problems  were  solved  with  the  codes  developed  here. 

The  initial  condition  for  both  problems  was  the  potential  flow  about  the 
body.  In  the  first  problem  we  assumed  that  the  boundary  layer  dev  'ops 
in  a  quasi-steady  way  but  the  vortex  sheets  develop  in  time  and  grow, 
feeding  their  influence  back  into  the  freestream  and  the  boundary  layer. 
In  this  formulation  we  allow  the  line  of  separation  to  relax  and 
relocate.  This  is  therefore  an  interactive  method.  In  the  second 
problem  we  specify  the  line  of  separation  either  arbitrarily,  or  we  use 
experimental  data;  for  example,  the  data  of  Meier  et  al.^^  No  boundary 
layer  solution  is  then  necessary.  However,  the  solution  is  again 
developed  iteratively,  by  allowing  the  vortex  sheet  to  grow  as  time 
increases. 

The  first  solution  was  constructed  as  follows.  A  set  of  skin 
friction  lines  based  on  the  exact  potential  flow  solution  was  calculated 


first.  For  a  prolate  spheroid  with  axes  ratio  1:4,  the  result  is  shown 
in  Fig.  9a.  The  strength  of  the  nascent  vortices  was  then  calculated 
and  allowed  to  drift  with  the  local  velocity.  The  inviscid  flow 
problem,  namely  Eqs.  (5)  were  solved  again  to  determine  the  new  loop 
circulations  on  the  body.  With  the  updated  potential  flow,  a  new  set  of 
skin  friction  lines  was  calculated  to  determine  a  new  envelope,  namely 
the  new  position  of  the  line  of  separation.  Once  again  the  stength  of 
the  nascent  vortices  was  calculated.  The  free  vortices  were  then 
convected  again  by  increments  proportional  to  the  local  velocities.  The 
process  was  then  repeated  for  a  total  of  six  time  steps. 

Figures  9b-e  display  the  development  of  the  vortex  sheet  after  the 
thrid  and  through  the  sixth  time  step  for  a  prolate  spheroid  at  an  angle 
of  attack  of  a  =  30“ ,  started  impulsively  from  rest.  In  these  plots, 
the  thick  lines  on  the  body  are  the  skin  friction  lines,  and  the  thin 
lines  are  the  inviscid  streamlines.  The  line  of  separation  is  clearly 
identified  in  all  figures  as  an  envelope  of  the  skin  friction  lines. 

The  separated  vortex  sheets  are  shown  in  the  form  of  a  grid  which  at  the 
nth  time  step  is  made  up  of  n  ribbons. 

It  is  clearly  desirable  to  allow  the  calculation  to  march  to  higher 
values  of  time.  It  is  believed  that  six  time  steps  are  not  enough  for 
the  full  development  of  the  wake.  However,  at  larger  time  steps 
numerical  instabilities  develop  which  destroy  the  smoothness  of  the 
wake.  Figure  10  displays  the  displacement  of  the  line  of  separation  as 
the  flow  interacts  with  the  boundary  layer.  It  is  seen  that  a  boundary 
layer  calculation  based  on  the  potential  outer  flow  predicts  very  early 
separation  in  the  forward  portion  of  the  body.  This  is  the  universal 
finding  of  researchers  who  employed  a  wide  variation  of  boundary- layer 


calculations  as  shown  in  Fig.  11.  It  should  be  noted  that  our  simple 
method  for  the  calculation  of  a  three-dimensional  layer  with  a 
potential,  unseparated  outer  flow  is  in  excellent  agreement  with  many 
other  complex  methods. 

In  Fig.  12  we  display  the  direct  comparison  of  flow  visualizations 
with  our  calculation.  The  streaklines  in  this  figure  were  obtained  by 
dyes  released  in  the  windward  side  of  the  model. The  dyes  are  aligned 
more  or  less  with  the  skin  friction  lines.  Moreover,  they  turn  sharply 
along  the  line  of  separation  and  lift  off  along  the  vortex  sheet.  It 
should  be  emphasized  here  that  the  first  vortex  line,  the  nascent  vor¬ 
tex,  is  located  at  a  distance  0.2R  above  the  body.  As  a  result,  viewed 
perspectively  from  an  angle  of  <>  =  90%  the  vortex  line  is  projected  on 
the  separation  line  only  if  its  polar  coordinate  is  i  =  90°.  For  larger 
ti's,  there  is  a  discrepancy  proportional  to  (0.2R)sin(ii,  as  shown  sche¬ 
matically  in  Fig.  13.  This  explains  partially  the  deviation  in  Fig.  12 
between  the  vortex  line  and  the  separation  line  at  the  leading  part  of 
the  body,  where  both  curves  tend  to  the  windward  side,  i.e.,  larger 
values  of  4. 

A  few  other  configurations  were  tested  with  our  interactive  code  as 
well.  Here  we  present  two  examples:  an  ellipsoid  with  axes  ratios 
1:6:4  and  a  projectile  shape.  Results  for  the  ellipsoid  at  angle  of 
attack  of  a  =  30°  are  displayed  in  Fig.  14,  representing  the  first 
instant  as  well  as  the  fifth  time  step  after  the  impulsive  start.  In 
Fig.  15,  we  show  the  results  of  the  interactive  code  for  the  flow  over 
an  ogive-cylinder  body  at  an  angle  of  attack  a  =  24.5°.  Figure  16 
displays  resuls  for  the  same  body  but  with  a  line  of  separation  fixed  at 
the  a  =  90°  meridional.  At  the  bottom  of  this  figure,  we  present  the 
calculated  vector  fields  at  two  cross-sectional  planes. 


To  facilitate  comparison  with  the  experimental  results  of  Meier  et 
al./°  it  was  necessary  to  prescribe  the  line  of  separation  because 
their  boundary  layers  were  turbulent  and  our  approximate  code  cannot 
handle  this  situation.  The  calculated  wake  is  shown  in  Fig.  17,  while 
in  Fig.  18  we  present  a  comparison  between  the  analytical  and  the 
experimental  velocity  fields  at  a  cross  section  of  the  flow.  Apparently 
the  vortex  sheet  has  not  developed  enough  after  six  time  steps.  The 
experimental  field  indicates  larger  velocities  in  the  vortex  as  well  as 
a  vortex  core  which  is  situated  a  little  closer  to  the  plane  of 
symmetry. 

A  tape  with  the  pressure  distributions  measured  by  Meier  and  his 
group  has  been  available  in  the  United  States  and  was  employed  by  other 
investigators  as  well  for  comparison.  In  Fig.  19  we  display  along  two 
meridionals  the  experimental  data  obtained  from  the  tape  and  the  present 
numerical  results.  In  the  same  figure,  the  potential  flow  results  are 
shown. 

Our  analysis  indicates  that  it  is  very  hard  to  predict  with  a 
discrete-vortex  method  the  pressure  distribution  undeneath  the  separated 
flow  region.  Apparently  the  wake  vortices  break  down  quickly  into 
turbulence  which  smooths  out  the  pressure  signature  on  the  body.  This 
became  more  obvious  when  we  plotted  our  predicted  distribution  as  a 
function  of  the  azimuthal  angle.  Another  reason  for  this  failure  could 
be  the  limitation  in  the  size  of  the  grid  and  the  number  of  vortex  lines 
which  make  up  the  wake.  An  example  of  pressure  distribution  at  x/L  = 
0.4812  is  shown  in  Fig.  20.  However,  this  relatively  good  agreement 
could  not  be  reproduced  at  all  stations  and  was  found  to  depend  upon  the 


size  of  the  time  step. 


It  is  believed  that  the  method  can  provide  reliable  results  for  the 
flow  and  pressure  distribution  over  the  attached  flow  region.  However, 


the  base  pressure,  which  incidentally  varies  with  the  length  of  the 
body,  can  also  be  predicted  simply  by  assuming  that  from  separation  and 
on  to  the  leeward  side,  the  pressure  is  uniform  and  equal  to  its  value 
at  separation.  Results  calculated  in  this  way  contrasted  with 
experimental  data  are  shown  in  Figs.  21  through  24. 

Conclusions 

The  present  calculations  have  further  demonstrated  the  potential  of 
the  vortex-lattice  method.  For  the  first  time,  a  boundary  layer 
solution,  although  somewhat  crude,  was  coupled  with  the  vortex-lattice 
solution  to  determine  the  location  of  separation  and  its  displacmenet  in 
time  as  the  wake  grows  and  rolls  up.  It  was  demonstrated  that  the 
strength  of  the  shed  vortex  can  be  determined  by  simple  approximate 
formulas  which  essentially  represent  integrals  of  boundary- layer 
profiles.  The  code  we  develooed  is  not  limited  to  the  configuration 
considered  here.  In  fact,  viscous  solutions  already  have  been  generated 
for  a  variety  of  ellipsoids  with  different  axes  ratios. 
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Velocity  distribution  along  the  meridional  of  a  sphere.  The 
vortex  lattice  had  6  circumferential  and  27  axial  panels. 
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Velocity  distribution  along  the  meridional  of  a  prolate 
spheroid  at  zero  angle  of  attack.  The  vortex  lattice  had  4 
circumferential  and  28  axial  panels. 
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g.  17  The  wake  over  a  prolate  spheroid  with  axes  ratio  1:4  and  the 
separation  line  fixed  according  to  Meier's  data. 
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Fig.  19  Experimental  and  analytical  pressure  distributions  for  the  case 
of  Fig.  17. 
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Fig.  21  Experimental  and  analytical  pressure  distribution  for  the  case 
of  Fig.  17,  along  the  axial  location  x/L  =  0.3501. 
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Fig.  22  Experimental  and  analytical  pressure  distribution  for  the  case 
of  Fig,  17,  along  the  axial  location  x/L  =  0.4812. 
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Fig.  23  Experimental  and  analytical  pressure  distribution  for  the  case 
of  Fig.  17,  along  the  axial  location  x/L  =  0.6061. 


Fig.  24  Experimental  and  analytical  pressure  distribution  for  the  case 
of  Fig.  17,  along  the  axial  location  x/L  =  0.7725. 


